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^ I The g-Gaussians are discussed from the point of view of variance mixtures of 

O ' normals and exchangeability. For each — oo < q < 3, there is a g-Gaussian distri- 

. bution that maximizes the Tsallis entropy under suitable constraints. This paper 

\^ , shows that g-Gaussian random variables can be represented as variance mixtures of 

normals. These variance mixtures of normals are the attractors in central limit the- 
' orems for sequences of exchangeable random variables; thereby, providing a possible 

PLh ■ model that has been extensively studied in probability theory. The formulation pro- 

r~| ! vided has the additional advantage of yielding process versions which are naturally 

I g-Brownian motions. Explicit mixing distributions for g-Gaussians should facilitate 

applications to areas such as option pricing. The model might provide insight into 
the study of superstatistics. 
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1. Introduction 

^ . Developments in nonextensive statistical mechanics based on an entropy 

[ proposed by Constantino Tsallis (1988) gave rise to the g-Gaussian distri- 

butions which are being applied in numerous research areas, see Gell-Mann 
and Tsallis (2004), Boon and Tsallis at al (2005), and TsaUis (2009). Ap- 
plications are based on experimental, computational, and analytical results. 
5^ I By definition, for —oo < q < 3, the g-Gaussian density takes the form 

gq{x) = Cq[l — (1 — q)x'^Y^^^^'^^ where Cq is specified in Section 2. Theo- 
retical under-pinnings of the g-Gaussians are founded on a novel g-algebra. 
The purpose of this paper is to propose a stochastic model within the frame- 
work of usual algebra that is consistent with many of the phenomena people 
are attempting to model by the g-Gaussian distributions with tails that decay 
as a power of x (those with 1 < g < 3) as well as the Gaussian distribution 
(g = 1). The model not only facilitates theoretical investigations, but helps 
to clarify why g-Gaussian distributions arise in certain physical applications. 

For each — oo < g < 3, Tsallis (1988) defined a g-Gaussian distribution to 
be that distribution which, under certain conditions, maximizes the Tsallis 
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g-entropy: 
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where e N is the total number of possible micropscopic configurations, 
{pi} are the associated configuration probabilities with Yl^iPi — 1' ^ 
a conventional positive constant. When q = 1 the definition is understood via 
the limit, in which case the 1-entropy recovers the Boltzmann-Gibbs entropy. 
The 1-Gaussian is simply the usual Gaussian distribution. The g'-Gaussian 
distributions have tails that decay on the order of when 1 < g < 3 

and have bounded support when g < 1. 

To a large extent, prevalence of the Gaussian distributions in theory and 
applications stems from their roles as attractors, e.g. via the classical central 
limit theorem for sequences of independent, identically distributed random 
variables with finite second moments. The dependence structure that should 
lead to g-Gaussian limits has remained elusive. The phrase "globally corre- 
lated" is often attached to phenomena in statistical mechanics to which Tsallis 
entropy is applied. 

There have been efforts to understand nonextensivity and the role of q- 

Gaussians from a probability point of view. For example, Umarov ct al (2008) 
proves g-central limit theorems using g-algebra and a notion of g-independence. 
Since the system developed is surprisingly nice mathematically, it is interest- 
ing to understand in the sense of usual probability the intuitive meanings of 
g-algebra and g-independence. Using usual algebra and usual notions of de- 
pendence, our paper establishes that g-Gaussians, with 1 < g < 3, have a 
role as natural attractors since they are a subset of the possible limits in a 
general central limit theorem for dependent random variables. This provides 
an explanation for the wide occurrence of at least these g-Gaussians. 

Our initial investigations were stimulated by the common features of two 
specific mathematical models for which it is possible to analytically verify 
whether or not the weak limit is g-Gaussian. Vignat and Plastino (2007) 
suggests a model which combines randomness in the normalizer with each 
i.i.d. random variable, and shows that the weak limit is g-Gaussian. Their 
model can be viewed differently, namely as a special case of a sequence of 
exchangeable random variables. Marsh et al (2006) provides practical models 
by applying Leibnitz triangles, in which the triangular arrays of random vari- 
ables are rowwise exchangeable, however weak limits are not g-Gaussian; see 
Hilhorst and Schehr (2007). Both examples are discussed further in Section 2 
since they motivated our consideration of exchangeability as a possible prob- 
ability model that might be connected to some applications being considered 
using Tsallis nonextensivity theory. 
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From a different but related direction, this paper is also stimulated by a 
possible connection between the proposed model and superstatistics. Beck 
and Cohen (2003) introduces the concept of superstatistics for generalized 
Boltzman factors derived from systems that evolve in complex environments. 
The expression for a superstatistic can be interpreted as a variance mixture of 
normals (VMON) (or scale mixture of normals) after standardization. More- 
over, in a specific case, the superstatistic reduces to the Tsallis statistic (Beck 
(2001)). VMONs are the attractors in a central limit theorem for exchangeable 
sequences, see Jiang and Hahn (2003). Furthermore, the reason behind study- 
ing superstatistics, such as a non-homogeneous background where a mechani- 
cal system expands, strongly suggests exchangeability in modeling the under- 
lying system. De Finetti's theorem ( Chow and Teicher (1978)) characterizes 
an infinite sequence of exchangeable random variables as a mixture of i.i.d. 
sequences, i.e. an infinite exchangeable sequence is conditionally an i.i.d. se- 
quence. Thus, if a mechanical system evolves in a complex, non-homogeneous 
background with the usual central limit theorem (for i.i.d. sequences) holding 
in each homogeneous part of the background, then the limiting process has 
marginals which are mixtures of normals. Further discussion of superstatistics 
appears at the end of Section 2. 

Similar arguments can be made in finance. Viewing market movements 
as responses to different specific information (interest rate, macro financial 
data, earnings, etc.), then conditioning on each specific type of information, 
yields a market which follows a random walk. Consequently, distributions that 
can model asset returns could naturally be mixtures of normals. Empirical 
evidence (Hall ct al (1989), Gribbin ct al (1992), and Kon (1984)) has shown 
that discrete VMONs fit asset returns data better than stable distributions 
or than the Student model. Stables have been studied more extensively than 
discrete VMONs. Moreover, research hterature on applications of continuous 
VMONs seems to be quite hmited. 

A drawback to applications of VMONs has always been that the mixing 
distributions are not generally easy to obtain from data. Expectation Max- 
imizing (EM) algorithms are usually employed. However, when the number 
of Gaussians in the mixture is large, or in many cases infinite, computations 
become very complicated. This paper shows that g-Gaussians are VMONs 
when 1 < g < 3 (and not when g < 1). Furthermore, the mixing distributions 
are calculated explicitly; thereby, facilitating application of these models. For 
example, the option price obtained from a VMON is the mixture of the option 
prices when the variance is fixed. Therefore, if one believes that g-Gaussians 
are good models for financial returns, then their option prices can be calculated 
simply by mixing the prices obtained from the classical Black-Scholes model 
with the known mixing distributions. Furthermore, studying mixing distribu- 
tions might provide insight into the "background" . For instance, a behef that 
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information drives stock prices suggests viewing the mixing distributions as 
models for the impact of information. 

The examples discussed above are each special cases of the model proposed 
in this paper. Our model of g-Gaussian distributions based on exchangeability, 
which is restricted to 1 < g < 3, is specified using usual algebra, rather than 
g-algebra, and has the following features: 

1) specification of g-Gaussian random variables as specific variance mix- 
tures of normals; 

2) associated central limit theorems for dependent exchangeable random 
variables with the g-Gaussians as natural attractors; 

3) associated stochastic processes which in a natural sense are g-Brownian 
motions. 

Within this model, as the latter two features indicate, the notion of "global 
correlation" is specified as exchangeabihty. Prom this point onward we con- 
sider the appearance of q-Gaussians only as variance mixture of normals under 
the exchangeability concept, not discussing their appearance from other con- 
cepts. 

1.1. Organization of the paper. Section 2 makes the connection between 
g-Gaussians and variance mixtures of normals as well as obtains a complete 

specification of the mixing distributions, which facilitates applications. Su- 
perstatistics are discussed as an example. Section 3 provides a detailed math- 
ematical description of exchangeability and establishes central limit theorems 
for exchangeable sequences and rowwise exchangeable triangular arrays. The 
examples motivating our consideration of exchangeability are discussed further 
in this section. Our concept of g-Brownian motion is introduced in section 4 
followed by some comparisons with a stochastic process in Borland (1998) that 
has g-Gaussian marginal distributions. Section 5 is the conclusion. 

1.2. Conventions. Throughout the paper, Z will be used specifically for a 
standard normal random variable that is independent of any other random 
variables. , stand respectively for convergence in distribution, 
in probability, and almost surely. It is worth re-emphasizing that all further 
discussion relies on usual algebra rather than g-algebra.. 

2. g-GAUSSIANS AND VARIANCE MIXTURES OF NORMALS 

The q- Gaussian distributions form a one-parameter family of distributions 
for — oo < g < 3 with densities specified by gq{x) — Cq[l — (1 — g)x^]^/(^~*\ 
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where the normahzing constant Cq is given by 



C: 
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h Til + y')^^dy 



(3-g)VT^r(5^)' 



— oo < g < 1, 

q = i, 

l<q<3. 



When q — 1, the expression for gq is understood by taking hmits and yields the 
standard normal density function. The support of the g-Gaussian is (— oo, oo) 

1 1 



if 1 < g < 3 and the compact set 



, if — oo < g < 1. 



A variance or (scale) mixture of normal distributions (VMON) by definition 
has characteristic function of the form 0(f) = exp[—t'^u/2)dH{u), with H 
a distribution function on [0, oo). The corresponding density is 



(27rM)-i/2 exp(-a;V(2M))rfif(M) 



Each VMON has a representation of the form VZ where V > a.s., and V 
is independent of Z. The VMONs include many commonly used distributions 
such as the symmetric stable distributions, the Cauchy, Laplace, double ex- 
ponential, logistic, hyperbolic, and Student distributions and their mixtures, 
plus many others. Sec Keilson and Steutel (1974) and Gneiting (1997). 

The connection between the g-Gaussians and the VMON can be made by 
applying the theory of completely monotone functions. A function h is com- 
pletely monotone on (0, oo) if and only if {—iy'h'^'^\x) > for x > 0, and 
71 = 0,1,2,... 

Andrews and Mallows (1974) shows that a symmetric density function fx{x) 
is a variance mixture of normals if and only if fxi^/x) is completely monotone. 
It is easy to show that gq{y/x) is completely monotone for g > 1 and not 
completely monotone for g < 1, which leads to the following statement. 

Theorem 2.1. q-Gaussians are variance mixtures of normals when 1 < g < 3 
and not mixtures of normals when g < 1. 

Since not every VMON is a g-Gaussian, it remains to identify the mixing 
distributions that yield the g-Gaussians for 1 < g < 3. In Beck (2001) the 
g-Gaussians are identified as where has a distribution with g de- 
rived from the degrees of freedom. Below we provide a direct proof using the 
Laplace transform technique when the inverse Laplace transform of the den- 
sity is known. It illustrates a method applicable to identification of unknown 
mixing measures from other variance mixtures of normals as well. 



Theorem 2.2. The q-Gaussian density for 1 < q <?>, can he expressed as the 
following variance mixture of normal densities: 

9,{x) ^ fv.z, = 1 ^ exp fyiv) dv 

where the mixing measure 

dH{v) = fv{v) dv = Cy^gCxp ( — 7777 — tv~9 ) "^^^^ dv, 



2(q-l)v 



Proo/. Let L denote the Laplace transform. Since L (exp(—6t) • t") — ^ ~'~ ^ 



for a > —1, the choice of a = — 1 = ^, b = yields 

q — I 9 — 1 

V V 9-1/ / + + 

Equivalently, 

poo 

C;[l - (1 - = / ex.p{-xH)dH{t), 

Jo 

where dH{t) = exp (^-^"^^ " and = T (^^— ^ j - 1)^- The 

following routine calculations verify the claim where the following substitu- 
tions are used: — > u in the second equality, 1/v ^ u in the third equality, 

and ^1 (^x"^ + u ^ v^ in the fourth equality: 
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= Cv,g I J- exp ) exp ( -7^7—^7^ ) v 9-1 



x/2^v "V2t;V V 2(g-l>;2 

2 ^g^) 



exp(— ^-(-(a^^H -)] dv 



2 V 1 
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It is not hard to verify that 

which completes the proof. ■ 

In Theorem 2.2, y is 1 over the square root of a distribution with the 
number of degrees of freedom being ^ — 1. Also notice that when q — 2, the 
g-Gaussian is Cauchy and 

When q < 4/3, EX"^ < oo. This situation is important in finance or risk 

management where the variance is often used to quantify the risk. In that 
situation, the g-Gaussians are superior to non-Gaussian stable distributions 
which fail to have finite second moments. When q < 3/2, EX < cxd. This 
situation is important in finance or insurance when measuring mean returns 
or mean expenses. 

The mixing distributions arc a type of generalized inverse Gaussian distri- 
bution, and the g-Gaussians are a type of generalized hyperbolic distribution. 
It is worth noting that both generalized inverse Gaussian and generalized hy- 
perbolic distributions are infinitely divisible. Since infinitely divisible laws are 
the only attractors for triangular arrays of independent random variables, this 
fact might help explain why g-Gaussians should be often observed without 
self- normalization (i.e. without pre-processesing). 

2.1. Example: Superstatistics. As a test particle moves from cell to cell, 
its velocity v satisfies a Langevin equation 

V = -jv + aBt, (1) 

where Bt is a standard Brownian motion. Instead of 7 and a being determin- 
istic, as in classical statistical physics. Beck and Cohen (2003) let /? = 7/(7^ be 
a random variable with density /(/3). In this setting, the widely used classical 
Boltzman factor e~^^ takes a generahzed form called a superstatistic: 

POO 

B{E) = / me-^'^dp, (2) 
^0 
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where E represents the energy of a microstate associated with each cell. Notice 
that the superstatistic in (2) is exactly a variance mixture of normals after 
standardization. The special case where k{f3) (the density of in our 

Theorem 2.2) is the density of a random variable yields the Tsallis statistic 



(Beck (2001)). 

3. Central limit theorems for exchangeable random variables 



A brief introduction to exchangeability is required. For details, see Chow 
and Teicher (1978). 

A sequence of random variables n > 1} defined on some probability 

space is said to be exchangeable if for each n, 

P{Xi < Xi, ...,Xn < Xn) = P{X^{1) < Xi, ...,X^(n) < Xn) 

for any permutation tt of {1, 2, n} and any Xi & H, i — 1, n. 

By de Finctti's Theorem, an infinite sequence of exchangeable random vari- 
ables is conditionally i.i.d. given the o"-field Q of permutable events. Further- 
more, there exists a regular conditional distribution for Xi given Q such 
that for each a; e Q the coordinate random variables {^n = > ^ > 1} (called 
mixands) are i.i.d. Hence, for each natural number n, any Borel function 
/ : R" R, and any Borel subset of R, 



The mixands are allowed to be on a different probability space than the X^s. 
For exchangeable sequences, the dependence never dies in contrast to weakly 

dependent sequences. Hence if the covariancc exists, it does not change along 
the sequence. Moreover, for an infinite exchangeable sequence the covariance 
is always non- negative. 

3.1. Example. The model proposed in Marsh et al (2006), considers iden- 
tical and distinguishable, but not necessarily independent binary subsystems. 
Let rjv,n be the probability that there are n subsystems in state 1, which is 
given by the Leibnitz rule: 



Since construction of the model only considers the number of subsystems in 
state 1, the order of I's and O's does not matter, a typical property of ex- 
changeable sequences. However, it is unclear whether each row of A^ 0-1 
random variables can be embedded in an infinite sequence of exchangeable 
random variables. Rodriguez et al (2008) and Hanel et al (2009) construct 



WITH Q- Gaussian limits 



P{f{x,, ...,x^)eB)= [ P{f{x,, ...,x^)e B\g)dP 




fN,n + fN,n+l — fN-l,n- 
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exchangeable models that can be embedded into an infinite sequence of ex- 
changeable random variables, again using the Leibnitz rule. 

Turning to central limit theorems, the following simplified and adapted 
version of Theorem 2.1 in Jiang and Hahn (2003) suffices for the needs of this 
paper. The proof is provided for completeness. 

Theorem 3.1. Let n > 1} be an infinite sequence of exchangeable ran- 

dom variables where Xi has mean if the mean exists or is symmetric other- 
wise. Assume < E < oo a.s. Then either 

i) ^ ^^V,-Z 

or 

n 

where Vi = ^/Var{^l) = ^Var{Xi\g) and V2 = E = E{Xi\ G). 
Proof. By de Finetti's Theorem, for any real x, 

<^= f ( T:=.i^-nEi, ^ ^^^^ ^ \ 

where cr(^i) denotes the standard deviation of ^1. If E^i = a.s., then the 
classical central limit theorem holds for the mixands and case i) of the theorem 
holds. If E$^i is not almost surely 0, replace the \/n in the previous equation 
by n. Then the first summand inside the regular conditional probability on the 
right side integral converges to zero almost surely and case ii) of the theorem 
holds. ■ 

Remark 3.2. The condition put on the mixands < E < 00 a.s. is 
generally weaker than assuming Xi 's have finite second moments. 

Remark 3.3. Case i) of Theorem 3.1 shows that q-Gaussians when 1 < q < 3 
are among the possible limits in a central limit theorem for exchangeable se- 
quences. However, they are not the only ones. It seems more natural to 
consider q-Gaussians as limit distributions in case i), since theoretically the 
range of distributions in case ii) is vast. For example, let ei's be i.i.d. with 
mean and standard deviation 1, and Y be any random variable that is inde- 
pendent of all Ei 's. Then {Y -\-ei, i > 1} is a sequence of exchangeable random 
variables and J20^ + ^i)/?^ ^• 

Remcirk 3.4. In case i), EX1X2 — if it exists, which means that q- 
Gaussians can be attractors of uncorrelated but not independent exchangeable 
sequences of random variables. In the second case, EX1X2 if it exists and 
in this case the central limit theorem is really the strong law of large numbers 
for exchangeable random variables. This case explains why the examples in 
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Rodriguez et al (2008) and Hanel et al (2009) achieve q-Gaussian limits with 
normalizers at the rate of n. 

Remcirk 3.5. Theorem 4-1 of Vignat and Plastino (2007) says that ifYi are 
i.i.d. in the domain of normal attraction of Z .and has a 'distribution 
(correcting a typo), then Xli converges in distribution to a q-Gaussian 
distribution. Their theorem in one dimension is a special case of Theorem 3.1 
in this paper with Xi = Vi ■ Yi where Vi = K However, in general, a sequence 
of exchangeable random variables does not necessarily have the form V -Yi. 

The next theorem is a triangular array version of Theorem 3.1, which seems 
to be new. 

Theorem 3.6. Let {Xn,i,i < n < n,n — 1,2,...} be a triangular array 
of rowwise exchangeable random variables that can be embedded into infinite 
exchangeable sequences. Assume Xn/t 's are centered or symmetric with < 
^ ^n,i < in probability for each n, and Var{^n^i) — >• in probability when 
n — > OO; with Vi > a.s. Then 



when \/nE^n,i ~^ ^2, o-nd 



V^-Z + V2 



with a > 1/2, when n^~°'E$^n,i V3 as n — > 00. 
Proof. By de Finetti's Theorem, for any real x, 



< X 



^ . /EkK^^^(^ ) + < ^"j dP. 

Vn J Jn \ V^-o-(Cn,i) / 

Using the definition of weak convergence in probability (Hahn and Zhang 
(1998)), and the conditions that v^^6t,i V2, f7(6t,i) Vi with Vi > 
a.s., and that the mixands have finite second moments almost surely, yields, 
when n — > 00 

P <x^ ^ J P'^iViZ + V2< x)dP = P{ViZ + V2< x). 

When n^'^'E^n,! ^3 with a > 1/2 and 1^2 7^ a.s., again, by de Finetti's 
Theorem, 

Since 
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when a > 1/2 and cr(^„,i) V2, 



and the second part of the theorem is proved. ■ 

Remark 3.7. When Xn^i 's are identically distributed, then V2 = and a = 1, 
which is consistent with Theorem 3.1. 

4. g-BROWNiAN Motion 

The hterature on g-Gaussians thus far has failed to allow the construction 
of a g-analogue of Brownian motion. Our representation of g-Gaussians for 
g > 1 in terms of variance mixtures of normals, using usual probabilistic and 
algebraic notions, leads naturally to the definition of a process that might 
naturally be called g-Brownian motion. However, we first require some defi- 
nitions. 

Definition 4.1. A stochastic process is called exchangeable if it is continuous 
in probability with Xq = and such that the increments over disjoint intervals 
of equal length form an exchangeable sequence. 

Definition 4.2. X has conditionally independent, stationary increments, given 
some a-field, if both properties of the increments are conditionally valid for any 
finite collection of disjoint intervals of the same length. 

These two definitions are connected by a theorem of Biihlmann which char- 
acterizes exchangeable processes on [See e.g. Theorem 9.21 of O. Kallen- 
berg (2002).] 

Theorem 4.3. Let the process {Xt)t>Q he R^-valued and continuous in proba- 
bility with Xq = 0. Then X is exchangeable if and only if it has conditionally 
independent, stationary increments given some a-field. 

Using these definitions, we can now define a g-Brownian motion. 

Definition 4.4. A q-Brownian motion, Bj, with 1 < q < ?> is a stochastic 

process having the following properties: 

1) conditionally independent, stationary increments; 

2) all increments are q- Gaussian; 

3) a.s. continuous sample paths. 

Theorem 4.5. (Existence) For each 1 < g < 3, a q-Brownian motion exists. 

Proof. Let Bt be a standard Brownian motion. When g = 1, B] = Bt is 
a process satisfying all three conditions with the conditional cr-field being the 
trivial cr-field. For 1 < g < 3, let Vg be a random variable with the density 
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fq given in Theorem 2.2. Then = Vq ■ Bt has conditionally stationary 
and independent increments given the value of Vq. All increments are q- 
Gaussian. Furthermore, a mixture of processes with continuous sample paths, 
has continuous sample paths. ■ 

4.1. Example. Compcirison with the Borland process. Borland (1998) 
provides another process with g-Gaussian marginals. Let Yt be the log returns 
of stock prices that follow the stochastic differential equation dY^ = /idt+adflt, 
where fit evolves according to dQt — P{^t)^^~'^^^'^dBt. The evolution of the 
probability distribution P of Y^ is nonlinear according to the nonlinear Fokker- 
Planck (forward Kolmogorov) equation 



Solutions for P{x, t\ y) are given by g-Gaussians for each fixed t, as established 
in Borland (1998) or using the g- Fourier transform as in Umarov and Queiros 



Our g-Brownian motion has conditionally, stationary and independent in- 
crements, each of which is r/-Gaiissian. This is expressed in the fact that its 
transition density is expressed in the form 



where p{t, x\y,v) for each fixed v satisfies the linear Fokker-Planck equation 

d 1 ^d^p{t,x\y,v) 

d^^'^'\y^'^ = r — — ' 

and the occurrence of the transition probability for a particular v is weighted 
according to the density fv{v) specified in Section 2. 

The Borland process clearly differs from our g-Brownian motion. Even 
though it has stationary increments (Umarov and Queiros (2009), its incre- 
ments are not conditionally independent. It does have continuous paths and 
g-Gaussian marginals. It is important to further clarify the nature of the 
increments for the Borland process for financial and other applications. 



This paper provides a probabilistic model for g-Gaussian distributions with 
1 < g < 3 based on exchangeability as one of the possible notions of "global 
correlation." The model is consistent with many of the phenomena for which g- 
Gaussian distributions are being used based on empirical evidence. In particu- 
lar, g-Gaussian distributions are variance mixtures of normals when 1 < g < 3 
and not when g < 1. Explicit mixing distributions are provided, which should 
extend further application of g-Gaussian distributions. An explanation for 




(2009). 




5. Conclusion 
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the wide occurrence of these g-Gaussians is their role as attractors via central 
limit theorems for exchangeable sequences and triangular arrays. A natural 
g'-Brownian motion is defined which can be viewed as an alternative driving 
process to the Borland (1998) process. The increments of the two processes 
have different characteristics and thus model different phenomena. The paper 
also makes a connection with superstatistics which are variance mixtures of 
normals after normalization. The Langevin equations that yield superstatis- 
tics can be viewed as stochastic differential equations driven by exchangeable 
processes. 
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